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(57) A method is provided for processing multi-component, multi-offset electromagnetic data measured at at 
least one multi-component receiver (20), the data representative of electric and magnetic fields due to a 
source, the at least one multi-component receiver being disposed at a depth greater than that of the 
source. The method comprises decomposing the measured multi-offset electric and magnetic fields into 
upgoing and downgoing components (21); and formulating a noise removal operator (22) from the 
downgoing components and the properties of the medium surrounding the at least one receiver. The 
method is used in seabed logging in CSEM applications. 
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1 

Processing Electromagnetic Data 

The present invention relates to the processing of electromagnetic data. In particular, 
the present invention is concerned with the calculation of a noise removal operator that 
5 attenuates certain parts of an electromagnetic field. 

The electromagnetic seabed logging (EM-SBL) technique is a new hydrocarbon 
exploration tool based on electromagnetic data, and is disclosed in Eidesmo et al., 
(2002) "Sea Bed Logging, a new method for remote and direct identification of 

10 hydrocarbon filled layers in deepwater areas". The Leading Edge, 20, No. 3, 144-152 
and in EUingsrud et al., (2002) "Remote sensing of hydrocarbon layers by seabed 
logging SBL: Results fi-om a cruise offshore Angola", First Break, 21, No. 10, 972-982. 
EM-SBL is a special application of controUed-source electromagnetic (CSEM) 
sounding. CSEM sounding has been used successfully for a number of years to study 

1 5 ocean basins and active spreading centres. SBL is the first application of CSEM for 
remote and direct detection of hydrocarbons in marine environments. The two first 
successfiil SBL surveys published were offshore West Afiica (Eidesmo et al and 
EUingsrud et al above) and offshore mid-Norway, R0sten et al., (2003) "A Seabed 
Logging Calibration Survey over the Ormen Lange gas field", EAGE, 65 An. hitemat. 

20 Mtg., Eur. Assoc. Geosc. Eng., Extended Abstracts, P058. Both studies were carried 
out in deep water environments (greater than 1,000 metre water depth). 

The method uses a horizontal electrical dipole (HED) source that emits a low fi-equency 
electromagnetic signal into the underlying seabed and downwards into the underlying 

25 sediments. Electromagnetic energy is rapidly attenuated in the conductive subsurface 
sediments due to water-filled pores. In high-resistance layers such as hydrocarbon- 
filled sandstones and at a critical angle of incidence, the energy is guided along the 
layers and attenuated to a lesser extent. Energy refiracts back to the seabed and is 
detected by electromagnetic receivers positioned thereupon. When the source-receiver 

30 distance (i.e. the offset) is of the order of 2 to 5 times the depth of the reservoir, the 
refracted energy from the resistive layer will dominate over directly transmitted energy. 
The detection of this guided and refracted energy is the basis of EM-SBL. 



2 

The thickness of the hydrocarbon-filled reservoir should be at least 50m to ensure 
efficient guiding along the high-resistance layer 

The electromagnetic energy that is generated by the source is spread in all directions 
5 and the electromagnetic energy is rapidly attenuated in conductive subsea sediments. 
The distance to which the energy can penetrate into the subsurface is mainly determined 
by the strength and firequency of the initial signal, and by the conductivity of the 
underlying formation. Higher fi-equencies result in greater attenuation of the energy and 
hence a lower penetration depth. The fi-equencies adopted in EM-SBL are therefore 
10 very low, typically 0.25Hz. The electric permittivity can be neglected due to the very 
low frequencies, and the magnetic permeability is assumed to be that of a vacuum, i.e. a 
non-magnetic subsurface. 

In terms of numbers, a hydrocarbon-filled reservoir typically has a resistivity of a few 
15 tens of ohm-metres or more, whereas the resistivity of the over- and under-lying 
sediments is typically less than a few ohm-metres. The propagation speed is medium- 
dependent. In seawater, the speed is approximately 1,700 m/s (assuming a frequency of 
1 Hz and a resistivity of 0.3 ohm-m), whereas a typical propagation speed of the 
electromagnetic field in water-filled subset sediments is about 3,200 m/s, assuming the 
20 same frequency and resistivity of around 1 ohm-m. The electromagnetic field in a high- 
resistance hydrocarbon-filled layer propagates at a speed of around 22,000 m/s (50 
ohm-m resistivity and IHz frequency). The electromagnetic skin depths for these three 
cases are approximately 275m, 500m and 3,600m, respectively, 

25 The electromagnetic receivers may be placed individually on the seabed, each receiver 
measuring two orthogonal horizontal components and one vertical component of each 
of the electric and magnetic fields. The HED source consists of two electrodes 
approximately 200m apart, in electrical contact with the seawater. The source transmits 
a continuous and periodic altemating current signal, with a fundamental frequency in 

30 the range of 0.05-10 Hz. The peak-to-peak AC ranges from zero to several hundred 
amps. The height of the source relative to the seabed should be much less than the 
electromagnetic skin depth in seawater to ensure good coupling of the transmitted signal 
into the subsurface, e.g. around 50-1 00m. There are several ways of positioning the 
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receivers on the seabed. Usually, the receivers are placed in a straight line. Several 
such lines can be used in a survey and the lines can have any orientation with respect to 
each other. 

5 The environment and apparatus for acquiring EM-SBL data are illustrated in Figure 1 . 
A survey vessel 1 tows the electromagnetic source 2 along and perpendicular to the 
lines of receivers 3, and both in-line (transverse magnetic) and broad-hne (transverse 
electric) energy can be recorded by the receivers. The receivers on the seabed 4 record 
data continuously while the vessel tows the source at a speed of 1-2 knots. The EM- 

1 0 SBL data are densely sampled at the source side, typically sampled at 0.04s intervals. 
On the receiver side, typical receiver separation distance is approximately 200-2,000m. 
Standard processing and interpretation of the acquired data can be performed in the 
common receiver domain or in the common shot domain, as long as data are sampled 
according to the sampling theorem (see, for example, Antia (1991) '"Numerical methods 

15 for scientists and engineers", Tata McGraw-Hill Publ. Co. Limited, New Dehli). 

The EM-SBL data are acquired as a time series and then processed using a windowed 
discrete Fourier series analysis (see, for example, Jacobsen and Lyons (2003) "The 
Shding DFT", IEEE Signal Proc. Mag., 20, No. 2, 74-80) at the transmitted frequency, 
20 i.e. the fimdamental frequency or a harmonic thereof. After processing, the data can be 
displayed as magnitude versus offset (MVO) or phase versus offset (PVO) responses. 

The principal wave types in the EM-SBL survey are illustrated in Figure 2. The wave 
types of main interest for hydrocarbon mapping involve only a single reflection 12 and 
25 a single refraction 13 at the target. These are detected as upgoing events by the receiver 

« 

3. A problem that arises in electromagnetic marine surveying is that electromagnetic 
energy may travel from the source 2 to the receiver 3 along many paths. The direct 
wave 8 is a signal transmitted directly from the source 2 to the receiver 3. The direct 
wave dominates in amplitude at short source-receiver separations, but is strongly 
30 damped at larger offsets since sea water has a high conductivity. In shallow water, EM- 
SBL exploration is complicated by source-excited waves received at the receiver array 
as downward-traveling waves which have been refracted (wave 1 1 ) and totally reflected 
(wave 10) off the sea surface 5. The air wave 1 1 is the signal that propagates upwards 
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from the source to the sea surface, horizontally through the air, and back down through 
the water column to the receiver. Due to the extreme velocity contrast between water 
and air, the critical angle for total reflection between sea water and air occurs at almost 
normal incidence. For angles of incidence greater than the critical angle, total reflection 
5 takes place, and the air volume acts as a perfect mirror for upgoing energy. The surface 
reflection 10 has its geometrical reflection approximately mid-way between the source 
and the receiver. In terms of signal strength at the receiver, the sea surface boundary is 
an efficient reflector at small to moderate offsets and an efficient refractor at larger 
offsets. The waves traveling downwards interfere with the upgoing waves from the 
1 0 subsurface. 

Reflections and refractions from the sea surface represent a severe problem, particularly 
in shallow water electromagnetic exploration. If the sea surface reflections and 
refractions are not sufficiently attenuated, they will interfere and overlap with primary 
15 reflections and refractions from the subsurface. 

In general, the water layer introduces a number of additional unwanted events that may 
interfere and overlap with primary reflections and refractions from the subsurface. A 
noise removal operator for removing unwanted events will be described below. The 

20 noise removal operator may also be known as a designature and denoise operator and is 
effective at substantially attenuating or completely removing the effects of the water 
layer present above the plane of the receivers in a typical EM-SBL envirorunent. The 
operator is effective at removing from electromagnetic data all events associated with 
any interface above the level of the receivers or with any interface at the receiver level. 

25 The operator is also effective at attenuating or removing the effects of the source 
radiation from the data. 

All energy and events caused by the medium above the receiver level will be referred to 
as "noise". 

30 

In order to provide accurate information about the subsurface target, it is desirable to be 
able to identify and substantially attenuate the incident wavefield due to the source and 
the noise from reflected and refracted waves received at the receiver. An important part 
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of any method for attenuating the source and noise wavefields will involve 
decomposing electromagnetic energy acquired at the receiver into its upgoing and 
downgoing constituents. There are two known approaches for this, see Amundsen, L., 
2003, Method for Electromagnetic Wavefield Resolution (WO 03/100467), and co- 
5 pending British Patent AppUcation No. 0407696.4. 

United States Patent No. 4,168,484 discloses a method for determining continuous and 
discontinuous impedance transitions in various media. The method involves disposing 
a source of electromagnetic radiation vertically above a number of receivers. Signals 
10 due to the source and due to reflections of media interfaces are recorded at the receivers 
and used to compute the incident and reflected waves, the incident and reflected waves 
being deconvolved to obtain the reflection impulse response. The reflection impulse 
response can be integrated to give the impedance transitions. 

1 5 According to a first aspect of the invention, there is provided a method as defined in the 
appended claim 1 . 

Further aspects and embodiments of the invention are defined in the other appended 
claims. 

20 

It is thus possible to provide a method which permits substantial attenuation of source 
and other noise components in electromagnetic data analysis. 

For a better understanding of the present invention and in order to show how the same 
25 may be carried into effect, preferred embodiments of the invention will now be 
described, by way of example, with reference to the accompanying drawings in which: 

Figure 1 illustrates the environment and apparatus for the acquisition of EM-SBL data; 

30 Figures 2a and 2b illustrate types of wave present in a typical EM-SBL environment; 



Figures 3a to 3c fiirther illustrate the wave propagation present in a typical EM-SBL; 
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Figures 4a to 4c illustrate the geometry of the method of an embodiment of the present 
invention; 

Figure 5 is a flow diagram illustrating a method in accordance with an embodiment of 
5 the present invention; and 

Figure 6 is a block schematic diagram of an apparatus for performing the method of an 
embodiment of the present invention. 

10 Optimal processing, analysis and interpretation of the electromagnetic data recorded at 
the receivers during a typical EM-SBL survey ideally requires full information about 
the field. 

The electromagnetic field will obey Maxwell's equations. In order to solve Maxwell's 
15 equations, the behaviour of the electromagnetic field at material interfaces and 
boundaries in the earth must be specified. At material interfaces, the tangential electric 
and magnetic fields are continuous. Even though all three electric and three magnetic 
components may be recorded, it is sufficient to record the two tangential components of 
the electric field and the two tangential components of the magnetic field. The nomial 
20 components of the electromagnetic field can be determined from Maxwell's equations 
when the tangential components are measured and the surroimding media properties are 
known. 

Figure 3a illustrates a multi-component source and multi-component receiver 
25 electromagnetic survey. The source 2 is a horizontal electric dipole that transmits a 
low-fi^equency electromagnetic signal down through the underlying rock formations. 
Using such a source, it is in principle possible to perform a two-component source 
survey where two orthogonal experiments are generated separately: one with the dipole 
antenna in the inline direction and a second with the dipole antennae arranged in the 
30 cross line direction. For each experiment, multicomponent electric and magnetic field 
sensors on a plane or along a line record the electromagnetic field. The source 2 emits 
electromagnetic waves with an amplitude which depends on the direction of 
propagation. Likewise, the receivers 3 record the electromagnetic waves with a 
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sensitivity depending on the angle of incidence. The arrows and dots in Figures 3a to 3c 
indicate the orientation of the sources and receivers: in the horizontal plane and 
perpendicular to the plane, respectively. The first two wave diagrams of Figures 3a 
show a transverse magnetic source, and the third and fourth show a transverse electric 
5 source. Upgoing and downgoing waves are emitted from the source and the receivers 
measure both upgoing and downgoing waves without distinguishing. 

A method of processing acquired or artificially generated electromagnetic data is 
described below which enables cancellation of the overburden effect In 

10 electromagnetic recording such as EM-SBL, the overburden is the water layer above the 
receivers, including the seabed interface. The method described below requires no 
information about the medium above and below the receiver plane, except for the local 
electric permittivity, magnetic permeability and electric conductivity at the receiver. 
For EM-SBL data in particular, only information of the electric conductivity is required 

15 due to their low-frequency nature. 

The method follows from the electromagnetic reciprocity theorem which provides an 
integral equation relationship between two independent electromagnetic fields defined 
in a specified volume enclosed by a hypothetical or physical surface. The relationship 

20 between the two fields is governed by possible differences in medium parameters, 
possible differences in source distributions, and possible differences in boundary 
conditions. The reciprocity theorem gives an integral equation procedure for 
transforming fields recorded in the physical electromagnetic experiment with the 
overburden response present into fields that would have been recorded in the 

25 hypothetical electromagnetic experiment with the overburden response absent. 
Mathematically, this follows from the reciprocity theorem by choosing outgoing 
boundary conditions for the desired field on the receiver plane. 

The wave-equation method that eliminates the overburden response is described as 
30 Lorentz designature/denoise analysis. This method preserves primary amplitudes whilst 
eliminating all waves scattered from the overburden. It requires no knowledge of the 
medium below the reviever level or above the receiver level. In the case where the 
subsurface is anisotropic and horizontally layered, the Lorentz designature/denoise 
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scheme can be simplified and implemented as a detemiinistic multicomponent source, 
multicomponent receiver, multidimensional deconvolution of common shot gathers. 
When the subsurface is isotropic and horizontally layered, the Lorentz 
designature/denoise de-couples on the source side into transverse electric and transverse 
5 magnetic problems, where a scalar field formulation of the multidimensional 
deconvolution is sufficient. 

The method begins from the assumption that the source is located in a horizontal plane 
anywhere in the water column strictly above the receiver plane. Further, the receiver 

10 measurements must allow a field decomposition on the receiver side just below the 
seabed into upgoing and downgoing wave components. From the upgoing and 
downgoing waves at the receiver level, the reciprocity theorem is used to eliminate the 
water layer response. The recorded physical electromagnetic data can then be 
transfomied to the desired data that would have been recorded in a hypothetical 

15 electromagnetic experiment without the water layer. The source in this hypothetical 
experiment is chosen to be a point source of electric current with some desired 
signature. A magnetic source may also be chosen and is an extension of the present 
invention that the skilled person would know to undertake. This situation is illustrated 
in Figure 3b. Since the water layer is absent and the incident field due to the source is 

20 removed, there are no downgoing waves at the receiver. The effect of the physical 
source and its radiation characteristics have been removed. This may be considered as 
new data having been designatured and denoised by a multidimensional signature 
deconvolution process. 

25 The Lorentz designature/denoise method, using data decomposed just below the seabed, 
replaces the water layer with a homogeneous half space with properties equivalent to 
those of the seabed. The designatured/denoised data will not contain the incident field. 
This data is highly usefiil for further processing and interpretation. 

30 Alternatively, the data may be decomposed into upgoing and downgoing components 
just above the seabed. The situation is illustrated in Figure 3c. hi this case, the effect of 
the seabed is still present in the designatured/denoised data. The effect of the water 
column and sea surface have, however, been eliminated. Applying the Lorentz 
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designature/denoise scheme just above the seabed is less preferable then applying it 
below the seabed because reflections and refractions from the incident field due to the 
point source will be present in the modified data. If application of the decomposition 
just above the seabed is the only possibility, a possible solution is to follow the 
5 designature/denoise processing with a further up-down field decomposition below the 
seabed. 

The notation used in the remainder of the specification is set out below in Table 1. Bold 
face type is used to distinguish matrices and vectors firom their components. The 
10 summation convention for repeated indices is used. Repeated Latin subscripts range 
over the values 1, 2 and 3 whilst repeated Greek subscripts take the values 1 and 2. The 
Kroenecker delta fimction is used 
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as is the Levi-Civita tensor, with components 



e ijk = 0, if any of ijk are equal 



20 



otherwise 



ei23= €312= €231 =-e213--e321 -"€132- 1. 



Table 1 



c 



system matrix, 

electric-magnetic field vector, 
complex velocity, c"^ = jus = -co'^tj^, 
unit vector along the - direction 



E = iE^^E^,E^) electric field, 

= (£, , ) = *E^"^ + ^^^^ horizontal electric field components. 
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F 
G 



upgoing components of horizontal electric field, 
downgoing components of horizontal electric field, 

4x1 source vector, 
2x2 Green's tensor, 

2x2 Green's tensor for the special case when source and 
receiver depths are close, and lateral source coordinates 

are zero, Xs = 0, 
magnetic field, 

horizontal magnetic field. 



K 

L,L-' 
n 



upgoing components of horizontal magnetic field, 
downgoing components of horizontal magenetic field, 

volume density of electric current, 
volume density of magnetic current, 
4x4 composition/decomposition matrix, 

2x2 submatrix of Z, 

unit vector normal to surface, 



p = (p^,P2) = fc/iD ,p = \p\ horizontal slowness vector, radial slowness 



vertical slowness. 



R 

X — , X2 > ^3 ) 

Six) 



S.. 



^ijk 

K = {k^^k^) = (op 



CO 



q^yjC -Pi -pI, qi =4^'^ -pI 5^2 =V^"' - 

reflectivity of subsurface 
wave vector, 

variables of Cartesian coordinate system, 
Cartesian horizontal coordinates, 
Dirac delta function, 
Kroenecker delta function, 

Levi-Civita tensor (the altering tensor), 
horizontal wavenumbers, 

circular frequency. 
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electric conductivity, 
magnetic permeability, 
electric permittivity, 

complex electric permittivity, 2^ = ^(1 + 

transverse admittance per length of the medium, 
7 = cr - ia>e = -icoe, 

longitudinal impedance per length of the medium. 



spatial derivative; d- = 



Gradient operator. 



The wavenumber, which characterizes the interaction of the EM field with the physical 
properties of the medium and firequency, can be written as 

A: = X"^. + iK_ , 



where 



r 



1 + 



cr 



2 \ 



Q) S 



1/2 



±1 



/ 



1/2 



The imaginary part of the wavenumber leads to the attenuation of a propagating EM 
10 wave in space. The wavenumber can also be expressed as: 
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with complex permittivity £ defined by 



6^ = £: 1 + 



coe ) 
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so as to absorb the conductivity as its imaginary part. This allows a unified treatment of 
an EM wavefield in both conductive (cr ^ O) and non-conductive (cr = O) media. For 

very high frequencies, (o»a I e ^ the wavenumber is real and given as 



and its dependence on the electric conductivity is negligible. Conduction currents are 
much smaller than displacement currents and can be neglected. In this circumstance the 
EM field propagates as a wave without significant attenuation. The scalar Green's 
10 function associated with the EM field, obeying (v^ + k^)j = -^7tS{x - x') , has the 
well-known form 



G = — exp(/M), 
R 



where /f =U- jcl. 



1 5 For very low fi-equencies, co « a/e, as in the EM-SBL experiment, and the field is 
said to be diffusive. The squared wavenumber is pwely imaginary. 



20 and its dependence on electric permittivity is negligible. Displacement currents are 
much smaller than conduction currents and can be neglected. Setting i = (l + V^* 
the wavenumber is written as; 



A«(l + /)c, 



with real component 
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V 
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In this circumstance, the scalar Green's function associated with the EM field is 



G = — exp(//c/j)exp(- /cR) 
R 



5 Since k is real the wave varies sinusoidally and is attenuated with distance. In one 
wavelength, the attenuation of the field is Itt, 

For EM-SBL wavefield decomposition, the complex electric permittivity is independent 
of the electric permittivity, but depends on the electric conductivity as 



(O 



The magnetic permeability is set to that of free-space (// = //o = 4;r • 10 ^ H/m) 
which is representative of a non-magnetic water layer and seabed. The complex 
1 5 velocity is then 



c = 



20 



The phase velocity is given by c^f^^col Re(^), yielding 



25 



Conductivity, measured in Siemens per metre, (or its reciprocal, resistivity) of sea water 

depends on salinity and temperature and typically is in the range cr- I-SS/tw. The 

salinity varies from sea to sea, but most major oceans have 3.5 percent weight. At zero 
degrees Celcius, the resistivity is approximately 0.34 Q m, and the conductivity is 2.94 
S/m. Under these conditions and at a frequency of 1/4 Hz the phase velocity in sea 
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water is c^^^ « 922 m/s. The skin depth S , where the EM wave will be reduced in 
amplitude by a factor of l/e , is 



At a frequency of Va Hz the skin depth in the sea water example is J « 586 m. 



Defining the geometry for the integral equation 

10 A volume V may be defined by the closed surface 5 = Z + with outward-pointing 
normal vector n , as illustrated in Figure 4a. E is a horizontal plane surface located at 
depth z~ infinitesimally above the multi-component receivers located at depth level . 
The Cartesian coordinate is denoted by x = {x,Xj), where X^i^x^^i)' notational 
convenience, x^=z , The z-axis, which is positive downwards, is parallel to it . The 

15 ATj , -axes are in the E plane. To simplify the analysis it is assumed that the medium 
is homogeneous and isotropic at depth z^ and in a infinitesimal region below. The 
overburden is the region for which z<z^ and the subsurface is that for which z> z^. 
Both may be arbitrarily inhomogeneous and anisotropic. Sj^ is a hemisphere of radius 
R. 

20 

In an EM-SBL survey, recording takes place on the sea bed. Due to continuity of the 
horizontal components of the EM field across the sea bed, the receivers may be assumed 
to be just below the sea bed. In this case, S coincides with the sea bed, and the 
overburden is the water layer, including the sea bed. Further below, the case in which 
25 the receivers sit just above the sea bed will be considered. 

An integral relationship between the multi-component source and the multi-component 
receiver data in the physical EM experiment will now be derived, containing the 
scattering response of the water layer above the receivers and the desired multi- 
30 component source, multi-component receiver data with that scattering response 
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attenuated. The physical source is assumed to separately generate two orthogonal 
electric currents along the horizontal axes of the Cartesian coordinate system. The 
desired muUi-component data are those data that would be recorded in a hypothetical 
muhi-component EM experiment from two orthogonally oriented sources of electric 

5 current acting separately with equal signatures when the medium above the receivers is 
homogeneous, extending upwards to infinity, with parameters equal to those at the 
receiver depth level (i.e. the sea bed). Magnetic point sources may also be used, but are 
not discussed further here. The overburden is therefore an isotropic halfspace. The 
geology below the receiver level is the same in the physical and hypothetical EM 

10 experiments. 

The physical EM experiment has a configuration as illustrated in Figure 4a. The 
recorded //'* component of the electric field vector at receiver location , just below 
Z , due to a source oriented in direction v at center coordinate with unknown source 
15 strength and radiation pattern is denoted by . Likewise, the component of the 
magnetic vector is denoted by . The source and field variables for the physical EM 
experiment, denoted as "state P", are listed in Table 2 below. 

The desired wavefields, and , that it is proposed to solve for are the responses 

20 of the medium fi-om two orthogonally oriented sources of electric current with desired 
signature or wavelet a corresponding to the dipole moment when the medium above 
the receiver level is a halfspace with properties equal to those of the sea bed as 
illustrated in Figure 4b. Z is a non-physical boundary. The desired electric and 
magnetic vector responses are recorded at location x^ just below E for the point 

25 sources located at jc; on S . The source and field variables for this hypothetical EM 
experiment denoted as ''state H" are listed in Table 2 below. 

To establish the integral relationship between the physical state P and hypothetical state 
H, the hypothetical "state H" is introduced, with wavefields E^^ and H^^ being the 
30 reciprocal wavefields to the ones in state H, obeying the reciprocity relation 
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5 Thus, E and are responses at location x~ on surface S due to a point source of 

electric current, with signature a, oriented in direction ju at location just below E 
as illustrated in Figure 4c. Surface S is, in the desired state H, an artificial, non- 
physical boundary. The source and field variables for state H are listed in Table 2 
below. 

10 

Table 2 





State P 


State H 


State H 


Electric Current 




aS{x -x~)e^ 


aS{x-x^)e^ 


Magnetic current 


0 


0 


0 


Electric field 








Magnetic field 









15 

Reciprocitv Theorem 

Reciprocity is an important property of wavefields. The reciprocity principle for 
elastostatic fields was derived by Betti and extended by Rayleigh to acoustic fields. In 

20 EM wave theory, reciprocity was introduced by Lorentz. The electromagnetic 
reciprocity theorem gives an integral equation relationship between two independent 
electromagnetic wavefields defined in a volume V enclosed by a surface iS. The 
relationship between the two wavefields is govemed by possible differences in medium 
parameters, possible differences in source distributions, and possible differences in 

25 external boundary conditions on S . 

Maxwell's equations for electromagnetic wave motion in an inhomogeneous medium 
can be expressed as: 



1 



J 
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15 



V X H{x, (o) - rj{x, co)e{x, co) - j{x, co\ 

V X E{x, 6}) + ^(jc, o))H{x, o)) = iSr(jc, 6)) 



5 In a domain or volume V enclosed by the surface S with outward pointing normal 
vector n , two non-identical electromagnetic fields denoted by the fields for "state A" 
and "state B", respectively may be defined. The boimdary conditions for the fields are 
not yet specified. State A is defined as 



10 VxH"" -Tj^E"" =/^ 

and state B is given as 



It is well known that by inserting special vectors, here denoted by Q , into Gauss' 
20 theorem. 



different Greenes vector theorems that are usefiil for studying wave propagation 
25 problems can be obtained. For EM waves, the specific choice 



Q^E^ xH^ -E^ xH 



30 is useful. Applying standard rules of vector calculus to V-g, yields the simple 
expression 

V.e = J/^(Vx£'^)-£:^ ^(VxH^)-H^ •(VxE^)+E^ .(Vxi/^) 
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Inserting this into Gauss' theorem leads to 

5 

^^dSn\E' xH' -E^xH']= \ydV[K^ -K' ^J^ -J^ E"" 
-[^^ -C'')h^^H' ^{tj' -rj')E'' 'E']. [1] 

Equation 1 is Green's vector theorem. It is also known as the reciprocity theorem, or 
10 integral representation, or integral equation for EM waves. The reciprocity theorem 
gives the relationship between two vector wavefield variables which characterize two 
states that could occur in the same domain or volume V. Each of the states may be 
associated with its own medium parameters and its own distribution of sources. On the 
right-hand side of Equation 1, the four first terms represent the action of possible 
15 sources in V. The two last terms under the volume integral represent possible 
differences in the EM properties of the media present in the two states. On the left-hand 
side of Equation 1, the surface integral takes into account possible differences in 
external boundary conditions. 

20 Reciprocity between state P and state H 

/\ 

The physical (state P) and hypothetical (state H ) EM experiments are described above 
and depicted in Figures 4a and 4c with volume V and enclosing surface 5 = 2 + 5^. 
Refering to the discussion of the previous section, state A is identified with state P 

A 

25 (Figure 4a), and state B with state H (Figure 4c). In both states, S is a plane surface 
infmitesimally above the receiver plane, and is a hemisphere of radius R . The field 
variables and sources for these states are defined in Table 2 above; thus giving in 
volume V for state A = P : 



30 E'=E^{x,cv\ 



f 



( 
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K* =0, 
5 /'*=0. 

The source term is zero since the soxirce, assumed to be a source of electric current 
oriented in direction v at center location , is outside V . Further, identify state 

B = H , so that in volume V : 

10 

7* =r]{x,a)\ 
15 =0, 

J* (jc) = a5{x — x^)e^. 

The fields are generated from a point source of electric current oriented in direction // , 
located at position x^ infinitesimally below surface Z . Inserting the above expressions 
20 into the reciprocity theorem yields 

Letting the radius R go to infinity, the surface Sg ->oo gives zero contribution to the 
25 surface integral. This is the Silver-MuUer radiation condition. Furthermore, taking into 
account that the surface Z is horizontally plane such that = -Sf^ and using that 

(CxZ)),=e^, CjD, 



30 gives 



20 



10 



Using the properties of the Levi-Civita tensor e^^ , gives 



=-l^ds{H,^E,, -H^,E,^ +^,^^2v -E^,H,^) 



Introducing the magnetic components = - and ^ = 1 in the above equation, 
the summation convention readily applies. This then gives: 



Equation 2 is the starting point for deriving the Lorentz designature/denoise scheme and 
describes the relationship between state P and state H and can be simplified by 
15 identifying proper boundary conditions for the fields on Z . hi the physical state P, 

and are siuns of upgoing and downgoing waves: 



= eV + E'^^K [3] 



20 

The physical fields, or equivalently, their upgoing and downgoing components, contain 
all information on the water layer overburden, including the effect of all physical 

sources. On the other hand, the data in the hypothetical state H experiment consist of 
upgoing events only, scattered fi-om the subsurface below S . In addition, the direct 
25 wavemodes fi*om the sources to the receivers are upgoing events since the sources are 

below the receivers. Thus in the hypothetical state HyE^^ and :7f„^ are purely upgoing 
fields: 



E = ^(^^ ; E^""^ = 0, [5] 




[6] 



Mathematically, to require outgoing (upgoing) boundary conditions on 2 for the fields 

of state H is equivalent to require the medium above S to be homogeneous. The 
boundary conditions of Equation 3 to 6 are most conveniently introduced into Equation 
2 by analysing the problem in the horizontal wavenumber domain, where upgoing and 
downgoing waves and their relation to electric and magnetic field vectors are 
analytically known. 

Relationships in the wavenumber domain 

A homogeneous isotropic region of the earth is now considered. Maxwell's equations 
can be written as a system of first-order ordinary differential equations of the form 



where the EM field vector B is a 4x1 column vector 

and the electric %^{Ey,E^y and magnetic H={-H^,hJ^ field vectors are 2x1 
column vectors. The 4x4 system matrix A is partitioned into four 2x2 submatrices of 
which the diagonal ones are zero, 



d^B=^icaAB-^F, 



0 



A 



A = 



A, 



0 



The symmetric submatrices A^^ and A2 , 



* 
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-1 



^1 

Pi Pi 



"PxPi 

^2 







PlPl 




.P\P2 
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are functions of the parameters in Maxwell's equations and of horizontal slowness 
When the source of magnetic current is zero {K -O), the source vector F is 



. 2. 



where 



F, = 



e P.J:, 


9 ^ 2 








J 2. 



For notational convenience, the explicit dependence of different quantities on 
frequency, wavenumber, depth, etc., is omitted. For instance, the electric field vector 
(E(jf,X3,co;jCs) recorded at depth due to a point source at location is in the 
wavenumber domain denoted £ or ^(Xj) with the understanding 
15 = ^e(x3 ) = fE(K:,X3 ,0); ) . 

Both the electric and magnetic field consist of waves travelling upwards {Jj) and waves 
travelling downwards (Z)). The electric and magnetic fields can then be expressed as: 



and 



25 The field vector B is decomposed into upgoing and downgoing waves of the electric 
field as 



I 



* 
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by the linear tranformation 



B = LfV, 



[7] 



10 



where L is the local eigenvector matrix of A (i.e., each column of L is an 
eigenvector). Equation 7 describes composition of the wavefield B from its upgoing 
and downgoing constituents. Given the inverse eigenvector matrix Z"' , the up-and 
donwgoing waves can be computed by evaluating 



This describes decomposition of the wavefield B into upgoing and downgoing waves 
15 of the electric field. 



The composition matrix 



I I 



[8] 



20 with inverse, the decomposition matrix. 



-1 



-L 



-1 

1 J 



25 



can be derived, where / is the 2x2 identity matrix, and 



1 

m 



/ _2 



P1P2 



\PlP2 9l J 



From Equation 7 and 8 it can be established that 



1 



* 
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From iV = L^B and the decomposition matrix L , the upgoing and downgoing 
electric-field components can be written as 



c(D)=l(«E_jr-i ^. 
2 



1 0 Similarly for the magnetic field: 



2 

1 5 In component form, the downgoing constituents are: 



£<1 



20 



1 



^9 



' 2 



1 



m 
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The corresponding upgoing constituents are: 



25 
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In a source-free homogeneous isotropic medium upgoing and downgoing waves satisfy 
the difTerential equations 

10 ^3^^"^ =-io)q^^. 



Making use of Parsevals' identity. Equation 2 yields: 



Introducing vector notation instead of using the summation convention, this can be 
20 written as: 



- (k, z; U.^)^,{-k,z;\x^\, 



25 where JC H2) = (-H^,Hx) and <E^ =={e,,E2) are the wavenumber domain 

magnetic field vector and electric field vector, respectively, and the superscript ^ 

denotes transpose. As detailed above, since the hypothetical state H fields and !E 
consist of upgoing wave modes only, they are related as 



26 



A (k) = jCj (k)e(k) , 



where i, , defined above as 



1 



^ 



<l2 P\Pl 



is a 2x2 matrix depending on the local medium parameters along the receiver spread. 
The matrix L, obeys the symmetry relation 



10 



£,(ic) = £[(«:)=L.(-/r). 



and fE^°^ are upgoing and downgoing horizontal components of the electric 
field £ respectively, such that 



15 



The physical state P fields DC and £ then are related as 



20 



[11] 



Inserting Equations 9, 10 and 11 into Parseval's identity, the upgoing waves £ 
cancel, so that 



(u) 



25 



[12] 



where the 2x2 matrix 



g-'=2L, 



« » 
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is interpreted as the inverse of the Green's tensor in a homogeneous medium when the 
source and receiver depths are infinitesimally close. Inverting this gives the Green's 



tensor 



1 r-I 1 

^ 2 leq 



r 



^1 

\rP\Pi 



P\P. 



J 



Furthermore, the vector 



10 



^(D)^[E(o),E{-)r 
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contains the elements of the downgoing wavemodes on each of the electric components 
and i?2 ■ Generally, for every shot location, fE^^^ may be calculated in the slowness 
(or wavenumber) domain from the electric and magnetic field vectors according to the 
downgoing components provided above and repeated here for convenience: 



' 2 



eq 



20 



^ 2 



1 



E2-Tr(piP2H2+qlH,) 
sq 



The scalars in front of the electric and magnetic field components are called 
decomposition scalars. The upgoing constituents are 



.(D) 



25 
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Eliminating the incident wavefield off the hypothetical state 



The desired field E^^ of the hypothetical experiment can be spHt into an incident wave 



field E^'^''^ propagating upwards fiom the source to the receiver, and the wavefield E^'^^ 
scattered upwards fi^om the subsurface. 



p — Mine) A(*c) 



In vector notation, the incident wave field, which propagates in a homogeneous 
medium, is the wavelet a multiplied by the Green's tensor G , that is. 



It can be further shown that: 



= afE<^^(- K, z^jjCs )exp(- iK - jr \ 



In Equation 12, on the left hand side the electric field can be split into upgoing and 
downgoing constituents and on the right hand side the hypothetical state electric field 
can be split into incident and scattered components. By identifying 



E^li^ k V. ) = ^ £ dK exp(- iK-Xr^^^\-K,t, \x, \ 



it can be seen that the downgoing part of the electric field cancels firom the left side of 
Equation 12, yielding 
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Using the reciprocity relation gives 



Using the property of the Green's tensor, G^"\k,z,\z~)=G^'^^{k,z;\z^), implying that 
£<™)(k,z, \z;)= (e("=^(k,z; \zJ, gives 



which can be interpreted as the "reflectivity" of the subsurface in the absence of any 
overburden. Given as linear combinations of £:J;7^ and E^^^K the elements of the 

reflection response are 




Finally, using Parseval's identity. Equation 13 reads in the space domain 

^rkK)= ldS{xyA^r\x)El^K^hl [14] 



where is the inverse Fourier transform of R^^, Equation 14 gives the sought-after 

integral relationship between the scattered field E^^^ (included in r^) in the 
hypothetical state H experiment and the state P total upgoing and downgoing fields 
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and . Thus, from the upgoing and downgoing wavefields, the reciprocity 

theorem has provided the theoretical basis for eUminating the physical response of the 
medium above the receiver plane (water layer overburden) in the multi-component 
source, multi-component receiver EM experiment. 

5 

Other than the position of the orthogonally oriented source elements, no source 
characteristics are required to eliminate all EM waves scattered from the overburden. 
Whatever the physical source characteristic is, it will be cancelled when solving for 

E^^'^ (or A"^) as this characteristic is present both at the left and right sides of Equation 

10 14 through the upgoing and downgoing fields. The multi-component sources have been 
transformed into point sources of electric current with the same frequency content as 
that of the physical source. This wave-equation method to eliminate the physical source 
radiation characteristics and waves scattered from the water layer overburden is denoted 
by Lorentz designature/denoise as the reciprocity theorem is originally credited to 

1 5 Lorentz. 

Equation 14 is a Fredholm integral equation of the first kind for the desired scattered 
fields, leading to a system of equations that can be solved for by keeping the 

receiver coordinate fixed while varying the source coordinate. Equation 14 can be 
20 compactly written as a matrix equation: 

^^^^ is found from the reflectivity r by multiplying in the wavenumber domain the 
25 relectivity R by the incident wavefield: 



10 



15 
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Wavenumber domain solution 

Fourier transforming Equation 13 over source coordinates receiver coordinates 

Xr yields the Lorentz designature/denoise procedure 



The leads to a system of equations that can be solved for and E^^^ by keeping the 

wavenumber conjugate to the receiver coordinate fixed while varying the wavenumber 
conjugate to the source coordinate. The coupling between the positive wavenumbers in 
the downgoing overburden response field with negative wavenumbers in the desired 
field (and vice versa) reflects the autocorrelation process between the two fields. In 
matrix form, the Lorentz designature/denoise process can be written as: 



Lorentz Decon volution: Horizontally Layered ID Medium 

An example of application of this method to a horizontally layered ID medium, 
constituting an embodiment of the invention, will now be described. For a horizontally 
layered medium, the response is dependent only on the horizontal distance between the 





source 




arbitrary horizontal vector. The shift variance implies 

^ap (-*^r l-*^) = ''^^ Of r + ^5 " Jf » \Xs^A' EQ^atiou 1 4 therefore can be written as 



that 



Making use of a variant of Parsval's identity yields 



4 
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X exp[/K -iXr+X, )J 



Fourier transforming with respect to Xr interchanging integrals gives 

5 

X tAt r ^Jfr exp[//, • (if - )1 



The last integral may be recognized as the Dirac delta function d{K-K/). Performing 
10 the integration over wavenumbers, using the Dirac delta function property 

r dKF(KpiK-K^)= F{k^\ where F{k) is any continuous function of k, and 
renaming jc, by #f , gives 

EfJ (ic, \x, )=R^(k,2, \x, , z; (k, z- \x, )exp(/if • x, \ 

15 

This can be written in terms of matrices as 

(E<") (k, z, |jCs ) = 'K^, z, Us > Zr" )^''\^ k, )exp(i • ) 
hiserting the expression for the reflectivity R (below Equation 13) gives 

20 



Solving for C^^*"^ gives the 'Xorentz deconvolution formula" 



25 £(-)(K,zJjC3,z;)=5E("H«,ZrkI«^°^(«,Zr"ks)r^''"'=^(«,z.|x..z;) [15] 



20 
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Equation 15 states that the desired scattered field is found by generalized spectral 
division between the upgoing and downgoing parts of the electric field, weighted by the 
incident wavefield of the desired state. The reflectivity of the subsurface can be given 
in terms of upgoing and downgoing constituents of the electric field as 



The Lorentz deconvolution can be expressed in terms of magnetic vector fields instead 
of electric vector fields. Using the relationships between upgoing and downgoing 
1 0 magnetic and electric vector fields given above yields 

where is the incident magnetic field in the desired state. Since Jf^"^^ is a 

15 downgoing field, it is related to S^^^ by 

^(inc) ^£^^i^^)^ 



Likewise, since Jf^^^ is an upgoing wavefield, it is related to (E^"''^ by 
ID isotropic medium 



As a fiirther example, a horizontally layered EM isotropic medium is considered. The 
25 wavefield is assumed to propagate in the a:i,X3 -plane such that p^—O, From 

Maxwell's equations two vmcoupled systems are obtained: one for E^^H^ waves, 
corresponding to EM waves with TM-polarization, and one for E^^H^ waves, 
corresponding to EM waves with TE-polarization. For TM-polarization the downgoing 
and upgoing waves are computed as 

30 
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' 2{ ' s ^ 



The electric dipole source is oriented along the -axis, giving the incident wavefield 



EtMXs ) = 2f g exp(- iK Xs)=^ exp(- in • Xs ) 



10 



The scattered part of the desired electric field is obtained according to Equation 1 5 by 
deterministic spectral deconvolution between the upgoing and downgoing part of the 
field itself: 



The scattered part of the desired magnetic field is correspondingly 



15 



^<r^(K,z,|;f,,z;)=[//<r>(ic,zJ:cJ///(f'(ir,z;|xJ]^r)(ic|;fJ 



20 



where the relationship 5£^^^ = yields 



Ht^ = - E^r^ = ^ exp(- iK ■ X, ) 
q. 2 



Multiplication by the incident wavefield is a signature process where the desired electric 
dipole source with wavelet a acts in the direction. 



25 



For TE-polarizatoin the downgoing and upgoing waves are computed as 



1 
2 



H, 



, ~ ■'-'2 ^2 ' 
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The electric dipole source is oriented along the jCj -axis, giving the incident wavefield 

5 The scattered part of the desired electric field is obtained according to Equation 15 by 
determisnistic spectral deconvolution between the upgoing and downgoing part of the 
field itself: 

10 

The scattered part of the desired magnetic field is correspondingly 
where similarly to the previous derivation, 

15 

h\^c) ^ ^ = - ^ exp(- iK . Xs ) 

/i 2 

Multiplication by the incident wavefield is a signature process where the desired 
electrical dipole source with wavelet a acts in the -direction. 

20 

The integral equation (Equation 14) can be modified to give a scheme for 
designature/denoise for electric field reflection data over 2D laterally inhomogeneous 
media. For TM-polarization with an electric dipole source oriented along the jc^-axis, 

25 ErAx,\x,)= [dS{x)-u(xM)E\f\x\x^ 

where r^, is the inverse Fourier transform of iijj , becoming 
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For the magnetic field, the corresponding designature/denoise scheme is 



where is the inverse Fourier transform of iiji ? becoming 



^21 21 A 5=0 J ^ 21 • 



10 



For TE-polarization with an electric dipole source oriented along the X2 -axis. 



1 5 where is the inverse Fourier transform of R22 , becoming 



2q 



n -Tl pise) _ rJ7 0''^)|y r F^^"^^ 

^22 " '^22 ~ L 22 |Aj=OJ ^22 * 



For the magnetic field, the corresponding designature/denoise scheme is 



20 



where r," is the inverse Fourier transform of R"^ , becoming 



R" = \h^""^ y V H^''^ 



25 
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Wavefield Decomposition Just Above Sea Bed 

The Lorentz designature/denoise method described above replaces the medium from the 
receiver depth level and upwards with a homogeneous overburden. In the previous 
5 sections, the receiver depth level was defined to be just below the sea bed by using the 
continuity of the horizontal components of the EM field across the sea bed interface. In 
this case, Lorentz designature/denoise processing gives idealised data without any 
events caused by the water layer and sea bed, 

10 However, instead of decomposing EM data into upgoing and downgoing waves just 
below the sea bed, the EM data can be decomposed just above the sea bed. In this case, 
the surface Z must be located infinitesimally above the depth of the wavefield 
decomposition. It follows that the Lorentz designature/denoise scheme replaces the 
water column and sea surface by a homogeneous water layer halfspace. This is 

1 5 illustrated in Figure 3c. Although the effects of the water column and sea surface are 
eliminated, Lorentz designature/denose processing will not remove any effects related 
to the sea bed. A disadvantage of applying the Lorentz designature/denoise scheme just 
above the sea bed is that reflections and refractions from the incident wavefield due to 
the point source of electric current will be present in the Lorentz designature/denoise 

20 data. These reflections will interfere with reflections and refractions from high- 
resistivity layers in the subsurface and may render the interpretation difficult. The 
solution to eliminate the sea bed reflection is to follow Lorentz designature/denoise 
processing with a fiirther up/down wavefield decomposition step below the sea bed. 

25 Redatuming 

The designature/denoised field as described above has been derived for the desired point 
source of electric current located just above the receiver plane. In marine EM-SBL the 
source is located a distance - z, above the receivers. The desired data can be 

30 redatumed to simulate acquisition from the physical source depth. Since the desired 
data are an upgoing wavefield, the redatuming is effected by multiplying the upgoing 
wavefield by a phase shift operator 
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The reciprocity theorem provides the theoretical basis for ehminating the physical 
response of a medium above a receiver level where EM waves are measured in a multi- 
5 component source, multi-component receiver experiment. The reciprocity theorem 
gives a procedure for transforming wavefields recorded in the physical EM experiment 
with the water layer overburden response present into wavefields that would have been 
recorded in the hypothetical EM experiment with the water layer overburden response 
absent. The transform process is called Lorentz designature/denoise. Other than the 
1 0 position of the sources, no source characteristics are required to eliminate all EM waves 
scattered from the water layer overburden. The radiation characteristics of the physical 
multi-component source are eliminated by a multidimensional source designature 
operation in the transformation from the physical experiment into the hypothetical 
experiment. 

15 

The Lorentz designature/denoise method requires that the physical wavefield is properly 
decomposed into upgoing and downgoing waves. Further the method of the 
embodiments requires no knowledge of the medium below or above the receiver level; 
and requires information only of the local and physical parameters along the receiver 
20 spread. The method additionally preserves primary amplitudes while eliminating all 
waves scattered from the water layer overburden. 

The Lorentz designature/denoise method is set out in the flowchart of Figure 5. At step 
20, EM data is acquired at at least one receiver. The data is then decomposed (step 21) 

25 into upgoing and downgoing components. The multidimensional designature and 
denoise operator that eliminates the response of the water layer overburden is computed 
at step 22 from the downgoing constituents of the multi-component data measurements. 
An integral equation is formulated at step 23 using the upgoing constituents of the 
multi -component field recording together with the multidimensional operator computed 

30 at step 22, and the desired source wavelet 23 for the electric current. The integral 
equation is solved at step 25 to give designature EM components with all of the waves 
scattered in the physical water layer overburden removed. 



> 
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In the case when the medium is anisotropic and horizontally layered, the Lorentz 
designature/denoise scheme greatly simplifies, and is conveniently implemented as a 
deterministic multidimensional deconvolution of common shot gathers (or common 
receiver gathers when source array variations are negligible). When the medium is 
5 isotropic and horizontally layered, the Lorentz designature/denoise decouples on the 
source side into TE and TM problems, with scalar field designature/denoise 
(deconvolution) processes. 

The schematic diagram of Figure 6 illustrates a central processing unit (CPU) 33 
10 connected to a read-only memory (ROM) 30 and a random access memory (RAM) 32. 
The CPU is provided with data 34 from the receivers via an input/output mechanism 35. 
The CPU then performs the wavefield decomposition 36, computes the signal removal 
operator from the downgoing components, and formulates and solves (numerically or 
analytically) the integral equation to provide the designatured data 37 in accordance 
15 with the instructions provided by the program storage 31 (which may be part of the 
ROM 30). The program itself, or any of the input and/or outputs to the system may be 
provided or transmitted to/from a communication network 38, which may be, for 
example, the Internet. 

20 It will be appreciated by the skilled person that various modifications may be made to 
the above embodiments without departing from the scope of the present invention as 
deifried in the appended claims. 
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CLAIMS: 

1. A method of processing multi-component, multi-offset electromagnetic data 
measured at at least one multi-component receiver, the data representative of electric 

5 and magnetic fields due to a source, the at least one multi-component receiver being 
disposed at a depth greater than that of the source, the method comprising: 

decomposing the measured multi-offset electric and magnetic fields into 
upgoing and downgoing components; and 

formulating a noise removal operator from the downgoing components and the 
10 properties of the medium surrounding the at least one receiver. 

2. A method as claimed in claim 1, comprising the further step of applying the 
noise removal operator to the measured electric and magnetic fields to attenuate the 
electric and magnetic fields due to the media at a depth less than that of the at least one 

1 5 receiver. 

3. A method as claimed in claim 1, comprising the further step of applying the 
noise removal operator to the upgoing components to attenuate the electric and 
magnetic fields due to (i) the media at a depth less than that of the at least one receiver, 

20 and (ii) the source. 

4. A method as claimed in any preceding claim, wherein the noise removal 
operator is formulated using electromagnetic wave theory. 

25 5. A method as claimed in claim 4, wherein the noise removal operator is formed 
using the electromagnetic reciprocity theorem between a first state and a second state. 

6. A method as claimed in claim 5, wherein the first state is the physical 
environment and the second state is a hypothetical environment in which the at least one 
30 receiver is bounded above by a homogeneous medium. 



7. A method as claimed in claim 6, wherein the homogeneous medium is free 
space. 
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8. A method as claimed in any one of claims 4 to 7, wherein the noise removal 
operator performs a multidimensional signature deconvolution process. 

9. A method as claimed in any preceding claim, wherein the decomposition of the 
measured data into upgoing and downgoing components is performed immediately 
beneath a horizontal plane in which the at least one receiver is disposed. 

10. A method as claimed in any one of claims 1 to 8, wherein the decomposition of 
the measured data into upgoing and downgoing components is performed immediately 
above a horizontal plane in which the at least one receiver is disposed. 

11. A method as claimed in any preceding claim, wherein the electromagnetic data 
is electromagnetic sea bed logging data. 

12. A method as claimed in claim 11, comprising the further step of redatuming the 
electromagnetic data using a phase shift operator. 

13. A method as claimed in any preceding claim, wherein the source emits 
multicomponent electromagnetic energy. 

14. A method as claimed in any one of claims 1 to 12, wherein the source emits 
single-component electromagnetic energy. 

15. An apparatus for processing electromagnetic data, comprising: 

a source for generating electric and magnetic fields; 

at least one receiver disposed at a depth greater than that of the source for 
measuring electric and magnetic fields; 

means for decomposing the measured fields into upgoing and downgoing 
components; and 

means for formulating a noise removal operator from the downgoing 
components and the properties of the medium surrounding the at least one receiver. 



• 
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16. A program for controlling a computer to perform a method as claimed in any 
one of claims 1 to 14. 

17. A program as claimed in claim 16 stored in a storage medium. 

5 

18. Transmission of a program as claimed in claim 16 across a commvmication 
network. 

19. A computer programmed to perform a method as claimed in any one of claims 1 
10 to 14. 



15 




Application No: 
Claims searched: 



GB04 14373.1 
1-19 



^^T. The 

• Intent • 
% QfiBce I 



43 



Examiner: 
Date of search: 



o 



INVESTOR IN PBOPLE 



Simon Colcombe 
14 October 2004 



Patents Act 1977: Search Report under Section 17 



Documents considered to be relevant: 



Category 


Relevant 
to claims 


Identity of document and passage or figure of particular relevance 


A 


All 


GB2385923 A 






(STATOIL) 


A 


All 


GB2333364 A 






(GECO) 


A 


All 


GB2296567 A 






(GECO) 



Categories: 



X Document indicating lack of novelty or inventive 
step 

Y Document indicating lack of inventive step if 
combined with one or more other documents of 
same category. 

& Member of the same patent family 



A Document indicating technological background and/or state 
of the art. 

P Document published on or after the declared priority date 
but before the filing date of this invention. 

E Patent document published on or after, but with priority date 
earlier than, the filing date of this application. 



Field of Search: 

Search of GB, EP, WO & US patent documents classified in the following areas of the UKC^ : 



GIN 



Worldwide search of patent documents classified in the following areas of the IPC 



07 



GOIV 



The following online and other databases have been used in the preparation of this search report 
WPI, EPODOC, JAPIO 



An Executive Agency of the Department of Trade and Industry 



